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SUMMARY i } 

Equations have been derived, using vector analysis, which pre- 
dict the incident energy at any given point on an earth satellite. The 
derivation is based on a non- spinning satellite and extended to includ 
spinning vehicles. The satellite was assumed to be oriented in any 
one of three ways — toward the sun, toward the earth, and tangent to 
the flight path at perigee. The amount of energy received depends 
upon the intensity of the radiation and the view factor between the 
object and source. The view factor equations are applicable to any 
mathematically describable shape. To determine the length of ex- 
posure to the sun, simple relations have been developed which pre- 
dict the ingress and egress points of the earth's shadow. The 
solution to the equations was programmed in a general manner on 
an IBM 7090. The input data required consisted of parameters to 
specify the orbit and an equation to describe the surface of the 
satellite. As an illustration a cylindrical satellite with hemispherical 
ends was chosen. The view factor equations were derived and 
numerically integrated, and the incident energy was calculated for an 
arbitrarily selected orbit. The results indicated that spinning reduces 
the peak heat fluxes in all cases, and the minimum incident energy 
occurred when the vehicle was oriented toward the sun. 

The methods outlined in this report are applicable to satellites 
of other celestial bodies also. The only additional information 
required is a knowledge of the position of the sun and the orbital 
plane relative to the body. 
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INTRODUCTION 


One of the major problems encountered in the operation of 
manned and unmanned space vehicles is heating by thermal radiation. 
The net heat received by the satellite may be controlled by many 
methods, both active and passive. For example, the vehicle may be 
oriented to obtain the minimum incident energy, or it may be shielded 
and insulated to reduce the heating, or a combination of these methods 
may be used. This report will consider the thermal advantages to be 
gained by minimizing the incident energy through vehicle orientation. 
For a spinning satellite, the problem is simplified somewhat since 
the incident energy is uniformly distributed about the spin axis. This 
fact will be used to determine the irradiation of a spinning satellite 
from the analyses of the non-spinning vehicles. 

An object in orbit about the earth receives significant amounts 
of radiant energy from three sources - direct solar radiation, earth 
radiation, and albedo or earth- reflected solar energy. As long as the 
satellite is near the earth it will always receive appreciable amounts 
of terrestrial radiation. However, it may or may not be receiving 
heat from the remaining two sources depending on the relative position 
of the earth, the sun, and the satellite. The amount of heat received 
from each source depends on the intensity of the radiation and the 
view factor between the object and the source. Equations for the view 
factors will be derived for any mathematically describable shape 
traveling in a circular or elliptical orbit. These equations will in- 
clude three possible vehicle orientations - earth-oriented, sun- 
oriented, and tangent to flight path at perigee. To determine the 
length of exposure to the solar heat source, equations will be derived 
which give the points of ingress and egress of the earth's shadow. 

The use of these methods will be illustrated by sample calculations 
for a cylindrical satellite with hemispherical ends and a verification 
obtained for certain simple cases. 
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DEFINITION OF SYMBOLS 

semi-major axis of orbit ellipse 

semi-major axis of earth’s shadow ellipse 

total surface area of radiant heat source 

total surface area of body receiving radiant energy 

finite area increments of satellite surface 

average albedo of a planet 

semi-minor axis of earth’s shadow ellipse 

number of days after vernal equinox 

diameter of a cylindrical satellite 

differential surface area on the radiant heat source 

differential surface area on the body receiving radiant 
energy 

differential area increment of irradiated surface 
differential radiant heat flux from dAi to dA 2 
eccentricity of orbit ellipse 

radiation view factor between the earth and the satellite 
area, A 2 i 

radiation view factor for a flat plate parallel to a 
planet's local horizon 

radiation view factor between the satellite area A 2 i 
and earth-reflected solar energy 


DEFINITION OF SYMBOLS (Continued) 


F s 

F v 


g 


h 


IR 

It 


Il2 


J. 


k 


k 2 


k 4 


L 

F 0 

LST 

N 

N, 

n 2 


- radiation view factor between the satellite area A 2i and 

the sun 

= radiation view factor for a flat plate vertical to the 
local planet horizon 

3 general mathematical function defined as 
g(X 4 , Xg, X 6 ) = 0 

altitude of a satellite above a planet’s surface 
= inclination of orbital plane to earth's equatorial plane 

- intensity of ear th- reflected solar energy 
= intensity of planetary radiation 

= total energy radiated per unit time and area by the 
earth 

intensity of radiation between the earth and the satellite 

= un ^t vectors along the coordinate axes of a rectangular 
cartesian coordinate system 

= 2. 206 X 10 4 (KM) 2 (constant) 

- 1.580570 X 10 9 (KM) 4 (constant) 

= length of cylindrical part of satellite 

- right ascension of sun in plane of ecliptic 
“ local sidereal time 

- North 

= unit vector normal to surface of heat source 
= unit vector normal to satellite surface 
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DEFINITION OF SYMBOLS (Continued) 
period of orbit 

semi-latus rectum of orbit ellipse 

incident radiant energy from body 1 to body 2 

incident radiant energy from earth to satellite area 
A *i 

incident energy from a planet 

earth-reflected solar energy incident on area A 2 
solar energy incident on area A 2 
radius from center of earth to the satellite 
radius of the earth = 6378. 150 ± 0.070 (km) 

radius of perigee 

radius from center of earth to any point on the earth 1 s 
shadow ellipse 

radius vector from earth surface element to satellite 
surface element 

solar constant 443 Btu/ft 2 hr 2% 

unit vector from center of earth directed toward sun 
location 

scalar components of solar vector along the X x , X 2 , X 3 
coordinate axis, respectively 

black body equilibrium temperature of the earth 

West 

a coordinate of a rectangular cartesian system 
rectangular cartesian coordinate system 


DEFINITION OF SYMBOLS (Continued) 


vector operator "del" 

right ascension of sun in plane of equator 

spherical angle in launch triangle (FIG 6) 

spherical angle in perigee triangle (FIG 6) 

angle in the plane of orbit, measured clockwise from 
perigee to the projection of the sun into the orbital 
plane (FIG 9) 

declination of the sun 

angle between unit surface normal to radiating surface 
and vector, r 12 (FIG Z) 

angle between the unit normal to the surface of the object 
receiving radiation and the vector, r\ 2 (FIG Z) 

angles between unit normals to the earth and sun 

angle between unit normals to the satellite and the sun 
(FIG 4) 

latitude of launch site 

latitude of initial perigee point 

obliquity of the ecliptic - Z3° Z7 ! 8. Z n 

angular coordinate on cylindrical satellite (FIG 15) 

geocentric angular coordinate for position description 
on the earth's surface (FIG Z) 

geocentric angular coordinate for satellite location 
(FIG 3) 
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DEFINITION OF SYMBOLS (Continued) 

= angular measurement in the plane of orbit (FIG 12) 

= geocentric angular coordinate for position description 
on the earth’s surface (FIG 2) 

= geocentric angular coordinate for satellite location 
(FIG 3) 

= angle subtending a side of launch spherical triangle 
(FIG 6) 

= angle subtending a side of launch perigee spherical 
triangle (FIG 6) 

= spherical coordinate (FIG 17) 

= true anomaly (FIG 8) 

= angle between the normal to the plane of orbit and the 
sun vector (FIG 9) 

= constant 3. 1416 

= spherical coordinate (FIG 17) 

= radius of planet 

= Stefan-Boltzman constant -• 174 5cl0~ 8 Btu/hr ft 2 (R°) 4 
= spherical coordinate (FIG 16) 

= spherical coordinate (FIG 1 6) 

= right ascension of ascending node 
= regression rate of ascending node 
= side of launch spherical triangle (FIG 6) 
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DEFINITION OF SYMBOLS (Concluded) 


00 

X 

argument of perigee 

00 

M 

advance of perigee 



Coordinate Subscripts 

1,2,3 

3 

earth centered triad (FIG 2) 

4, 5, 6 

3 

satellite coordinate system (FIG 5) 

7,8,9 

S 

system given by first rotation in equatorial plane 

10,11, 12 

S 

orbital plane triad 

13, 14, 15 

= 

earth-oriented triad 

16, 17, 18 

= 

system given by first rotation to align with the sun 
(FIG 10) 

19, 20, 21 

3 

sun-oriented system 

22, 23, 24 

= 

positive hemisphere system 

E.O. 

3 

earth-oriented system 

P 

3 

perigee system 

S. O. 

3 

sun-oriented system 
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RADIATION HEAT TRANSFER 

In this section, equations will be developed which predict the 
thermal energy incident on the satellite. The analysis will be made 
for a non-spinning vehicle and then extended to include the spinning 
satellite. The irradiation will be assumed to be independent of wave- 
length and angle of incidence. It is convenient to use the celestial 
sphere concept which assumes that the earth is stationary and the 
heavenly bodies revolve about it. Two rectangular cartesian co- 
ordinate systems will be employed: one moving with the satellite, 

but fixed relative to it, and a geocentric equatorial system. By 
calculating the maximum heating that could be received by an earth 
satellite, it will be shown that only the earth and sun contribute 
significant amounts of thermal radiation. 

Consider two diffusely radiating surface elements dAj and dA 2 , 
as shown in FIG 1, which are separated by a non-absorbing medium. 
The radiation emitted from dA* and striking dA 2 is 


dqi- 


cos Yi cos y 2 


L 12 


htz| 


dAi cLA-2 


( 1 ) 


where Ii 2 is the radiation intensity. The angles Yi and Y 2 are 
the angles between the respective normals and a line connecting the 
two elements. The total radiation received per unit time by area A 2 
from Ai is 


3- * / I 

Ai Az 


~ cos y. COS . v _ z _ dAi dA2 


1 ^ 2 ! 


( 2 ) 


Or assuming that A! is a black body at temperature Tj , then 


Qi- 


S 


COS Vl COS Y2 
2 


A, A? 


rTz 


dAj dA 2 (3) 
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Figure 1. Radiation Heat Transfer Between TWo Surfaces 
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and the limits of integration are cos Yi >. 0, cos Y 2 — 0 (4) 

A. EARTH PLANETARY RADIATION 

If area dAi represents an elemental area on the surface of 
the earth and dA 2 is the satellite area as shown in FIG 2, then 


cos Yi 


r 12 



ri2 


and 


cos Y 2 



( 5 ) 


( 6 ) 


where the unit normal to the earth's surface is 


Ni = sin rp 


cos 0] (ij) + sin n j sin Gj (i 2 ) + cos rp (i 3 ) 


(7) 


and the elemental area is given by 

dAi = Re 2 si n rp dtp d©! (8) 

The position of the satellite relative to the earth is given by the 
vector R, as shown in FIG 3, and is defined in terms of the geo- 
centric angles r| 2 and 0 2 . If the magnitude of R is denoted as R 
then the vector v^ z from the earth's surface to the satellite is 

r\ 2 = (R sin p 2 cos 0 2 - Re sin hi cos 9 i ) (r) + 

(R sin p 2 sin 0 2 - Re sin hi sin 0j ) (i 2 ) + (9) 


(R cos r| 2 


Re cos r| 4 ) (i 3 ) 
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Figure 2. Radiation Heat Transfer Between the Earth and Satellite 
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Figure 3. Satellite Position Co-ordinates 
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The area dA 2 on the satellite will depend on the particular shape of 
the vehicle and will be defined later. 

The only vector as yet undetermined is the satellite unit normal. 
Let X 4 , X 5 , X 6 represent the coordinate axes of a rectangular 
cartesian coordinate system which moves with the satellite. This co- 
ordinate system is fixed relative to the satellite and chosen to give 
the simplest mathematical description of the satellite surfaces. 

The general equation of the satellite surface would be 

g(X4, X 5 , X f ) = 0 (jq) 

and its unit normal at any point is 



(ID 


with the sign chosen so that N 2 is always directed outward from the 
surface . 

The integration of equation (3) may be simplified by studying 
the functional dependence of the vectors on area A 2 . r^ is not a 
function of position on the satellite surface, since even for a 100- 
foot satellite in a low orbit (100 miles), the differences in rq 2 at 
extreme points on the satellite are negligible. Also by subdividing 
the satellite surface into increments (A 2i ) over which the 
direction of the unit normal is constant, the cos y 2 is made piece- 
wise independent of A 2 . The equation for the energy received by 
A 2 i from the earth is 

o-T , 4 f 

Ql ~* 2i= 1T" J dA *i 

A *i 


or 


Qi-*2j & T i* A 2 ^ F 2 i 


; 


cos Vt cos y 2 


12 


dA, 


( 12 ) 


(13) 
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where F 21 is the view factor defined as 

dA x (14) 

Additional information concerning the definition of the view factor is 
given in Appendix A. 

To apply equation (13) the surface temperature of the earth must 
be determined. For a given point on the earth the temperature varies 
with the time of day, season and atmospherical conditions, thus 
making it impossible to use the Stefan-B oltzman Law. However, it 
is possible to compute the total energy radiated by the earth by 
means of an energy balance, independent of these variables, and 
obtain a mean effective temperature. The thermal energy incident 
on the earth consists almost entirely of solar radiation. Since the 
sun is so far from the earth, it can be considered as a point source 
emitting radiation which impinges on the earth in parallel lines. 

The intensity of this solar energy is called the solar constant, S, 
and is based on the irradiation of a flat plate normal to the sun's 
rays at the earth's mean distance from the sun. The net energy 
absorbed by the earth is the difference between the incident solar 
radiation and the fraction that is reflected as defined by the albedo. 
Since the mean temperature of the earth and its atmosphere does not 
vary appreciably over extended periods, it may be concluded that 
this absorbed energy is in turn reradiated. Both the solar constant 
and average albedo, A, are reasonably well known for the earth 
so that this reradiated energy can be readily evaluated. If It is 
the total energy radiated per unit area and time, the energy balance 
is 


21 


= fe . I 


cos Yi 


r 12 


(1 - A)S7TR E 2 = 4 7TR E 2 It 


(15) 


This assumes that all the surface of the earth is radiating with the 
same intensity as a black body; therefore 


It = "T , 4 


(1 - A) 


4 


S 


( 16 ) 
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The most accurate value for the solar constant is that given by- 
Johnson (Reference 1), as 


S 


443 ± 


9 


Btu 
ft 2 hr 


and the average albedo is 0 . 40 (Reference 2). Substituting these 
values into equation (16) gives a mean effective temperature for the 
earth of 443. 5°R, and an intensity of 66. 36 Btu/hr-ft 2 . 

B. SOLAR RADIATION 

Since the irradiation of a flat plate normal to the sun’s 
rays is known (solar constant), it is easy to find the incident energy 
on any flat surface. The amount of heat received will be proportional 
to the projected area in the direction of the sun. For example, if 
the unit normal to the area A makes an angle y 2S with the sun's 
rays, as shown in FIG 4, then the radiation received from the sun 
would be 

Q s . = S cos y 2S A 2i (17) 

From the form of the above equation, and the definition of a view 
factor, as explained in Appendix A, it can be seen that the solar view 
factor is 

F s = cos y zs = S • N 2 (18) 

where 

s = Si (iT) + s 2 (i 2 ) + S 3 (it) (19) 

C. ALBEDO RADIATION 

The third source of energy irradiating the satellite is earth- 
reflected solar energy. Neglecting the atmosphere, the amount of 
solar energy incident on the earth per unit time and area would be 
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S cos Yis (20) 

where 

— ► — ► 

cos Yis = s • Ni (21) 

However, all of this energy is not absorbed; some of it is reflected 
back to space. The average fraction reflected is defined as the 
albedo and accounts for the actual reflected radiation regardless of 
the cause of its reflection. Therefore, the amount of the energy 
returned to space, per unit time and area, by the earth would be 


A S cos Yis 


( 22 ) 


If the energy is reflected diffusely according to Lambert's Cosine 
Law then the intensity in any direction is 


!r 


A 5 cos Yis 

TT 


(23) 


The amount of this energy received by the satellite can now be 
computed in a manner similar to direct earth radiation. The only 
difference is the intensity of the radiation. Such an approach gives 


Q Ri = S AA 2i Frj 
where 


(24) 


FRi = x / 


COS Yl3 cos Yl COS Y 2 


\r~7z\ 


dA, 


and the limits of integration are 


(25) 


cos Y 2 > 0 


cos Yi .> 0 


cos Yis 2 0 


( 26 ) 



19 


The limits of integration given by equations (4) and (26) deter- 
mine whether or not a surface is receiving radiation. Specifically, 
the limit cos > 0 means that a given element of area on the 
earth is in a position to emit energy to a particular satellite location. 
If the cos y^ is positive, this means that the element of area on 
the satellite is in a position to receive the radiation. The limit on 
cos y insures that the earth surface area element is irradiated 
by the sun and therefore should be considered when summing the 
reflected solar energy* 

D. OTHER CELESTIAL SOURCES 

To determine whether or not the radiation from a nearby 
celestial body should be considered, the maximum possible heating 
loads were calculated* At any given altitude above a planet this 
maximum incident energy would occur on a flat plate which is 
parallel to the local planet horizon. The view factor between the 
plate and the planet is determined by Smolak* s equation (Reference 3) 
and found to be 



(27) 


The intensity of the planetary radiation is determined in a manner 
similar to that previously described on Page 16. The albedo values 
are taken from Reference 2 and listed with the calculated intensitis 
in Table 1. The incident thermal energy per unit area on a flat plate 
parallel to the planet* s local horizon would be 



Table 2 gives the amount of incident energy for various 
altitudes above the moon and the neighboring planets. It can be seen 
that planetary radiation is negligible beyond 30, 000 miles from the 
surface of the bodies considered. Therefore, the only sources which 
contribute significant amounts of thermal radiation to an earth 
satellite are the sun, the earth, and albedo or earth-reflected solar 
energy. 
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TABLE 1 


Planetary Constants 


Celestial 

Body 

Earth 

Moon 

Venus 

Mars 

Btu/hr-ft^ 

66.36 

102.86 

50.78 

40.48 

A 

0.40 

0.07 

0.76 

0.15 


TABLE 2 

Planetary Radiation to a Flat Plate Parallel to the Planet's Local Horizon 


h 


q (Btu/hr-ft^ 
P 

! ) 


(Statute 

Miles) 

Earth 

Moon 

Venus 

Mars 

100 

63.0 

86.4 

48.3 

36.7 

500 

52.0 

48.5 

39.7 

26.2 

1000 

42.1 

28.0 

32.0 

18.8 

5000 

13.0 

3.17 

9.48 

3.48 

10000 

5.33 

1.00 

3.87 

1.19 

20000 

1.88 

0.274 

1.30 

0.357 

30000 

0.906 

. ... 

0.636 


40000 

0.538 


0.386 
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E. SPINNING SATELLITE 

Spinning the satellite causes the heating loads to be uniformly 
distributed about the spin axis. To calculate the incident energy the 
satellite is divided into segments which are bounded by planes per- 
pendicular to the spin axis. The total irradiation of each segment is 
then calculated for the non- spinning case and this value divided by the 
exposed surface area of the segment. This solution assumes a 
constant spin rate during any one revolution, which will usually be 
the case. However, variable rates can also be accounted for by 
weighting the averaging process. 

ORBITAL PARAMETERS 


The parameters which describe the satellite's orbit and motion 
will be given in this section. Spherical trigonometry will be used 
to relate the launching conditions to the resulting orbit, and the 
movement of the satellite in the orbital plane will be described in 
polar coordinates. Because the earth is not a perfect sphere, 
orbital parameters will not remain fixed but will experience per- 
turbations, which will be accounted for by Krause's equations 
(Reference 4) . 

As the vehicle travels in orbit it is assumed to be oriented in 
any one of three ways: toward the earth, toward the sun, or tangent 

to the flight path at perigee. In each case the vehicle will be assumed 
to be capable of maintaining its orientation and not be influenced by 
drag or electromagnetic forces. The primary interest here is in 
the angular relations and not in expressing the coordinates of a point 
in terms of the different coordinate systems; therefore, the trans- 
position of the origins will be neglected. The equations relating 
various unit vectors of the different coordinate systems are given 
in detail in Appendix B, based on the analysis of this section. 

The position of the sun relative to the geocentric coordinate 
system will be derived and the attitude of the earth's shadow. 
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cylinder 1 will be determined. Since the intersection of the orbital 
plane and the shadow-cylinder describes an ellipse, the ingress and 
egress points will be found by simultaneously solving the equations 
for the orbital and shadow ellipses. Such an approach assumes that 
the satellite is a point, which is a reasonable assumption for the 
purposes of angular location. 

The orbital plane is defined by the angle, i, which it makes 
with the equatorial plane of the earth and the right ascension of the 
ascending node, f 2 , as shown in FIG 5. The position of the orbit in 
its plane is given by the argument of perigee, go, measured from 
the ascending node to the perigee point. 

The initial track of the satellite, on the earth's surface, is 
illustrated in FIG 6, based on the assumption of a spherical earth 
and a launch trajectory confined to the orbital plane. 

To correlate the orbital and launch parameters, spherical 
trigonometry is used. This method requires the solution of two 
similar right spherical triangles. These triangles are formed by 
the earth's equator and the track of the satellite, with the third 
side being the meridian of longitude through either the launch site 
or perigee projection. The launch triangle is solved first since two 
of its parts are known. However, the known quantities here are the 
inclination of the orbit, i, and the side opposite which is the latitude 
of the launch site. This case is ambiguous, giving two possible 
solutions; thus, a sketch should be drawn to interpret the results 
properly. The spherical angle relations are 


sin 6j^ *= 

sin i sin goj^ 

(29) 

tan = 

cos i tan goj^ 

(30) 


i Actually the earth's shadow is a cone with an included angle of 
0. 53° and an apex located 746, 800 nautical miles from the earth's 
center (Reference 5). However, the cylindrical assumption is 
practical since even at a distance of 10, 000 nautical miles from the 
earth the diameter of the cylinder would only be 80 nautical miles 
more than the cone, or one percent of the earth's diameter. 


Descending Node 



Figure 5. Illustration of the Orbit Orientation Relative to the Earth 
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and 


sin 3 l = sin Xl esc oo^ (31) 

Using the above results, the parts of the spherical triangle with 
the meridian of perigee as a side may now be found. The known 
quantities are the angle of inclination, i, and the argument of perigee, 
oo, which is given by 

.oo = oo L + downrange cutoff (32) 

The remaining relations are calculated by the following equations 
which are similar to equations (29) through (31) above 

sin 6 * sin i sin oo (33) 

tan X = cos i tan oo (34) 


and 


sin 6 * sin X esc oj (35) 

p 

The date and time of launch appear in the determination of £2, 
since 

a * LST - A L (36) 

as can be seen from FIG 6. The angle, Al, is found from the solution 
of the launch triangle, and Reference 6 gives the local sidereal time, 
LST. 

The polar form of the equation of an ellipse which is 
a( 1 - e 2 ) 


R * 


1 + e cos v 


(37) 
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will be used to describe the position of the satellite in the plane of 
orbit. The angle, v, is called the true anomaly and is measured 
counterclockwise from perigee. The radius is measured from the 
perifocal point to the satellite. In the case of an earth satellite the 
perifocal point is chosen at the center of the earth and the radius, 

R, as shown in FIG 3. The eccentricity of the ellipse is e, and a is 
the semi-major axis length. 

Because of the oblateness of the earth, the perigee and nodal 
locations do not remain fixed but vary according to the equations 
given below (Reference 4). The location of the ascending node will 
change with time according to 


Q - £2 \ initial passage 

of the satellite 


( 


Number of days since \ £2 

initial passage of the J (38) 
satellite 


where 


n = - 


¥ < 24) cos 1 [ 3 i|f 

(t - 1 3in i i)j 


* "sM' *H 


radians 

days 


(39) 


and 


K z = 2.206 X 10 4 (KM)‘ 


K 4 = 1.580570 X 10 9 (KM) 4 


Rp (1 + e) 


(40) 

(41) 

(42) 


= (2.764 X 10-6) 


Rp (KM) 
1 - e 


hours 


(43) 
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The perturbation of the argument of perigee is 


« = 4r (24) 


o K 2 
3 —7 

( 1 - 4 sin 2 i) 

+ 10 73&I 

(\ + J e 2 ) 

w 

\ 2 / 

w 

V 4 / 


^1 - 5 sin 2 i + sin 4 i^ j - (cos i) 


• radians 

(O) 


day 


(44) 


and the position of perigee will be 


) 


to = to \ initial passage + / Number of days \ 

of the satellite I since initial passage! 

\ of the satellite / 


(45) 


The relationship between the satellite coordinate system (X 4 , 

X 5 , Xfc) and the earth-centered system (X 1( X 2l X 3 ) depends on the 
attitude of the vehicle. The derivation of the transformation equations 
between these two systems will be explained by a series of rotations. 
(FIG 5 is an aid to visualize the physical significance of these 
rotations.) First, rotate the Xj, X 2 , X 3 coordinate system about the 
X 3 axis through the angle, £2. This rotation places the X 7 axis on 
the line of nodes. 


The matrix relating the new triad to the earth-centered system 
is given below in terms of the direction cosines. 


cos £2 cos (7r/2 - £2) 

cos (tt/2 + £2) cos £2 

COS 7T / 2 COS 7T / 2 




' " 


COS 7T /2 


Xi 


COS 7r/2 


X 2 

= , 

cos 0 


x 3 




. 4 

. 


X 7 


X 8 


X, 


or simplifying 


( 46 ) 
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cos £2 

sin £2 

0 


Xi 


X 7 

- sin Q 

cos £2 

0 


X 2 

= 

X 8 

0 

0 

1 


x 3 


X, 


( 47 ) 

To place the base of the triad in the plane of the orbit, the X 7 , X 8> 

X^ system is rotated about the X 7 axis through the orbital inclination 
angle, i (FIG 7 ) 


The matrix relating this system to the previous one is 


1 0 0 

0 cos i sin i 

0 - sin i cos i 




r 1 

X 7 


0 

X 8 

II 

X11 

x? 


X12 

. 


. 


( 48 ) 


Further orientation in the plane of orbit is accounted for by a 
rotation about the X 12 axis. When the coordinate system is rotated 
through the angle, to, the Xjq axis will be pointing toward perigee 
and the Xjj axis will be tangent to the flight path of the satellite 
at the perigee point. This coincides with one of the desired satellite 
orientations; and if the axis of the orbiting triad, which is tangent 
to the flight path, is X 5 , then the matrix relating the satellite axis 


ystem to the X 10 , 

X n , 

X,2 triad 

is 





COS CO 

sin co 

0 


X 10 


X4 


- sin co 

COS CO 

0 


X u 

= 

x 5 


0 

0 

1 


x 12 


X 6 


Now the satellite vectors may be expressed in terms of the geocentric 
coordinates by successive application of the transformation equations. 


To orient the satellite so that the X5 axis is always pointing 
toward the center of the earth requires only a turning of the satellite 
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in the plane of its orbit. The angle through which the satellite must 
be rotated is measured relative to the perigee attitude of the vehicle 

7 r 

and, as can be seen from FIG 8, is simply — + v. 


If the earth-oriented triad is called Xj 3 , X^, X 15 as shown in 
FIG 8 , then the relationship between it and the perigee system is 
given by the matrix 


- sin v 

- cos v 
0 


cos v 
- sin v 


0 





r 


X 4 


X ,3 


x 5 

II 

X14 


X 6 


X„ 

- 

. 


_ - 


(50) 


The sun orientation is more difficult to express mathematically 
since the satellite must not only be turned in the plane of its orbit but 
must also be elevated out of the orbital plane. These rotations will 
be expressed in terms of the angles which the sun's unit vector 
makes with the perigee system as defined in FIG 9. The scalar 
components of the sun vector in the perigee system are 


sin £, cos r S • U 

cos £ = S • 4 

and 

■ ► ■ -> 

- sin £ sin F - S • 4 


(51) 

(52) 


(53) 


Since the vectors in the above equations are known in the geocentric 
system (equation 19 and Appendix B),the indicated vector operations 
may be completed and a scalar expression obtained for the right hand 
side of each equation. Using these results and the fact that 0 _< £ <_ 7 T, 

equations (51-53) may be solved for £ and I""* . 
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/ 

Normal to the ^ / 

Orbital Plane e^/ 



Figure 9. Position of Sun Relative to Satellite Co-ordinate System 
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/i - T7 

By rotating the perigee system through the angle — + 1 , as 

shown in FIG 10, the X 5 axis is aligned with the sun’s projection in 
the orbital plane (FIG 9). The relationship between the new triad 
(X^, X 17 , X 18 ) and the perigee system is 


- sin T 

- cos r 

0 


x 4 


Xi6 

cos r 

- sin F 

0 


x 5 

= 

X 1 7 

0 

0 

1 


X 6 


X 18 


(54) 

To complete the alignment with the sun, the satellite must now be 

7T 

elevated about the X ^ axis through the angle, — - so that the 

X 17 axis now points directly at the sun. The axes of the triad 
resulting from the last rotation are labeled X^, X 20 , X 21 . The 
transformation between this system and the previous one is given by 
the following matrix array: 


1 

0 

0 


0 

sin | 

- cos £ 



r 




Xi6 


X 19 


x 17 

11 

X 20 


X 18 


X„ 

- 

„ 


. 


(55) 


To aid in remembering the orientation associated with each 
arithmetical subscript, the important coordinates will hereafter be 
referred to by more informative subscripts as listed below. 


(pe rigee system) 
(earth- oriented) 


and 


X 4 , x s , x 6 = Xp, y p , Zp (56) 

Xi 3 , X M , X 15 = X E. O. ’ YE. O. ’ Z E. O. 

(57) 

Xi 9 , X 20 . X 21 = x s o. » YS.O. > z S.O. 


( sun-oriented) 


( 58 ) 
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The location of the sun may be found by using the celestial sphere 
concept and spherical trigonometry as shown in Reference 7. The 
projection of the path of the sun on the celestial sphere is called the 
ecliptic and makes an angle, e , with the celestial equator of 23. 45° 
or 23° 27’ 8 . 2 ,T . Using the vernal equinox as a reference point the 
right spherical triangle to be solved is formed by the ecliptic, the 
celestial equator, and the meridian circle passing through the sun's 
position. If the right ascension in the plane of the ecliptic, Lq, was 
known, then the equations could be solved. The relations from 
spherical trigonometry would be 


tan a - tan Lq cos <e 


(59) 


and 


sin Y = sin e sin Lq 


(60) 


where a is the right ascension of the sun in the plane of the equator 
and Y Is the declination of the sun. The expression for Lq will be 
found by use of the equation relating true anomaly to mean anomaly, 
since the mean motion of the sun is well known. This equation is 
(Reference 4) 


5 

v = M + 2 e sin M + -7 e 2 sin 2M + . . . (61) 

4 

where the eccentricity of the earth's orbit is 0.0167 (Reference 5). 
Since the sun completes an orbit every 365.25 days, the mean daily 
motion would be 


360 

365. 25 


0.9856 

day 


(62 


and the mean anomaly expressed relative to the vernal equinox would 
be 
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M - 0.9856 (77 + Number of days after v. e.) (63) 

The number 77 appears because the passage of the vernal equinox 
occurs 77 days after perihelion of the earth. However, basing the 
position of the sun on the calendar introduces an error in the 
calculations. The calendar neglects the one-quarter day and com- 
pensates by adding a full day every four years. To locate the true 
position of the sun in any given year a correction factor is added to 
account for this discrepancy. The right ascension in the plane of the 
ecliptic would then be 


L ,0 = 0.9856 (D + 77) + 1.9481 sin ^0.9856 (D + 77)J 

+ 0 . 0207 sin^O. 9856)(D + 77) J - 


(64) 


77 617° + Number of years since leap year 

4 


where 77. 617° is the true anomaly of vernal equinox; and the number 
of years since leap year will be 0, 1, 2, or 3. A new leap year is 
not counted as four years since the last one, but as zero years. Now 
Lo can be calculated, and equations 59 and 60 solved for a and 7 . 
The relationship between the coordinates a and 7 , and the earth 
triad, is shown in FIG 11 along with the conversion equations. 

The intersection of the earth's shadow cylinder with the orbital 
plane is illustrated in FIG 12 and will be described in terms of 
previously defined parameters. The angle of intersection is measured 
by the angle between the normal to the orbital plane and the sun 

vector. The plane view shows elliptical shape of the locus of inter- 
section and its relation to the X 16 , Xi 7 , X 18 triad system. The 
equation of the earth's shadow ellipse is 



1 


(65) 




Earth-Sun Line 


Figure 12. Earth's Shadow 
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where 


a s “ 


RE 


cos 


and b s = Re (FIG 12). Or, using the angle 


9 and the radius R s , measured from the center of the earth to any 
point on the shadow ellipse, the equation may be expressed as 

Rs = 7 ~ ~ RE --Z t\ T — 7TT7 ” 


COS 


0(1- cos 2 £,) + cos 2 £ 


The position of the shadow ellipse in the orbital plane is defined 
relative to perigee by the angle F . 

The equations of the two ellipses in the orbital plane (the shadow 
ellipse and the satellite orbital ellipse) are now known. Positions on 
satellite orbit are specified by the true anomaly, v, while the angle 
0 locates points on the shadow ellipse. From FIG 1 3 it can be seen 
that these angular coordinates are related by the following equation 


v = 0 + tt/2 - F 


(67) 


The ingress and egress points are two of the four possible inter- 
sections of these ellipses, and may be found by simultaneously 
solving their equations. This has been accomplished on an electronic 
computer by simultaneously computing the two radii and comparing 
them. By limiting the comparison to the values of v corresponding 
to a 0 range of 0 <_ 0 <_ 7T , the proper intersections are obtained 
(FIG 14). 


EXAMPLE PROBLEM 


As an example, the incident energy on a cylindrical satellite 
with hemispherical bulkheads will be calculated for all three vehicle 
orientations. The thermal heating depends on the intensity of the 
radiation and the view factor between the object and the source. The 
intensities are known and the view factor equations may be evaluated 
by double integration. The input data needed to perform the integration 
consists of an expression for the satellite surface normal and specifi- 
cation of the orbital parameters. 




Satellite 

Orbit 



Figure 14. Ingress and Egress of the Earth's Shadow 
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The surface normal equations will be derived in detail for the 
satellite oriented tangent to the flight path at perigee. The other 
derivations are similar, and the results are given in Appendix C. 

Consider the cylinder oriented as shown in FIG 15. The equation 
for its surface would be 


(x p ) 2 + (z p ) 2 


d 2 

4 


( 68 ) 


and the unit normal would be (see equation 11) 


N 2 


(Ip) ( x p) + (kp)(zp) 

d/2 


(69) 


or using the coordinate transformation 

Xp = d/2 cos £ (70) 


Zp = d/2 sin t, 

The surface normal then becomes 

N 2 = (ip) cos £ + (kp) sin £ (72) 

The relationship between the satellite coordinate system and the 
earth-centered triad is given by equations (B4-B6) of Appendix B. 
Substituting those expressions into equation 72 gives 
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N, 


(7) 

cos C, 

|cOS GO 

cos f2 

- 

sin co 

cos i sin 

q \ 

+ 


sin £, 

^ sin i 

1 

sin ^2 

i) 

+ 




(S) 1 

|cos 

J cos OJ 

sin f2 

4 

sin co 

cos i cos 

°] 

+ 


- sin 

C [ sin 

i cos 


}♦ 




( 13 ) | 

|cos £ 

J sin ou 

sin i J 


+ sin 

£ |cos ij 

) 



(73) 


The hemispherical tank ends will be described by placing limits 
on the equation of a sphere. The end nearest the origin of the 
satellite coordinate system (negative hemisphere) would be part of 

a sphere of radius — whose equation is 


(xp) + (y p ) + (zp) 2 = d 2 / 4 


(74) 


and the unit normal would be given as (see equation 11) 


n* = (ip)( x p) + (Jp)(yp) + (kp)(z p ) 

d/2 


(75) 


To restrict the surface to a hemisphere, it is convenient to switch 
to spherical polar coordinates as shown in FIG 16. Then the unit 
normal is 


N 2 - sin 4 1 sin (J> (ip) + cos <f> (jp) + sin <f> cos ^ (kp) (76) 
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Figure 16. Co-ordinate System for Description of Negative Hemispherical 
Bulkhead 
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and the limits which define the hemisphere are 

0 <_ Ztt (77) 

*72 <_ <f> <_ 7r (78) 

To describe the other hemispherical end (positive hemisphere) 
the origin of the coordinate axes is shifted to x p = 0, y p = L, 

z p = 0, where L, is the length of the cylindrical portion of the 

satellite. Denoting this new system as shown in FIG 17, the unit 
normal here is 


N 2 - sin p sin p (i 22 ) + cos p (i^) + sin p cos p (i^) (79) 


and 




*22 


jp “ *23 



(80) 


The angular limits are 


0 < p < 27r 


(81) 


0 < p < tt/2 


(82) 


The equations for the hemisphere unit normals are expressed in 
terms of the earth-centered coordinate system just as the unit 
normal to the cylinder was giving 

N 2 - (ij) | s in 4* sin <|> Jcos oj cos £2 - sin o> cos i sin £2 J 

cos <(> sin u> cos £2 - cos gj cos i sin £2j 

cos ^ sin <j> £sin i sin £2jJ + 
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(i. 2 ) |sin 4 sin <J> Jcos co sin £2 + sin co cos i cos £2 J + 

cos <j> Jcos co cos i cos £2 - sin <o sin £2 j 

cos 4 sin 4 J- sin i cos £ 2 j | + ( 8 : 


( 13 ) | s in 4* sin 4 | sin <0 s in i J + cos 4 | cos w s in i j 
cos 4 sin 4 ^ cos i J j 


= (ii) |si 


in ^ sin p cos go cos £2 - sin <0 cos i sin £2 I + 


cos p - sin co cos £2 - cos co cos i sin £2 + 


cos H- sin p I sin i sin £ 2 j j + 

(ij) |sin H- sin p |cos co sin £2 + sin co cos i cos £2 j + 

r 1 < 84 > 

cos p cos 00 cos 1 cos 12 - sin oo sin J 2 + 

cos p. sin p | - sin i cos £2 j J + 

(13) |sin (x sin p |sin oo sin i j + cos p |cos oo sin ij + 

cos p sin p | cos ij j 

The orbit selected has the following pertinent characteristics 


i = 28° 45' 


e = 0 


a - 4118.7 miles 
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Launching will be assumed to be made from Cape Canaveral on 
May 11, 1962, in the south-east direction. The geographic co- 
ordinates of the Cape are 

latitude, 28° 45' N 
longitude, 80° 34' W 


From this information, the "launch triangle" may be solved using 
equations 29 through 31; the results are 


= 7t/2 


X L = Jr/ 2 


3 L - tt / 2 


(87) 


Assuming a downrange cutoff of 15° gives a value for the argument of 
perigee of 


oj = 105° 


from equation 32. The other parts of the perigee triangle are found 
from equations 33 through 35 and are 


6 = 

27° 41' 

(89) 

X = 

107° 

(90) 

3 = 

81° 55' 

(91) 


The local sidereal time at midnight on May 11, 1962, is given in 
Reference 5 as 


LST 


228° 31' 22" 


(92) 
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and the right ascension of the ascending node is then 
« = 228° 31' 22" - 90° = 138° 31' 22" 


(93) 


from equation 36. 

Since this is a circular orbit, equation 37 reduces to 


R = a 


(94) 


The perturbations of the orbital parameters are calculated from 

equations 39 through 44; these give the position of the ascending node 
as 


138 31' 22" + (-7.6725) (Numbers of days since 

initial passage of the (95) 

satellite) 


and the argument of perigee varies according to 


- 105° + ( 12. 423)(Number of days since initial ( 96 ) 

passage of the satellite) 


The ingress and egress points of the earth's shadow are found to 
beat v - 270.2° and 57.7°, respectively, for the initial orbit. 

The position of the satellite in terms of earth-centered coordinates, 
02 and r| 2 , is needed to compute the vector r^ 2 (equation 9). The 
equations for these angles were given with FIG 3 and are 


COS T|2 


2k 

R 


(97) 


sin 0 2 


X 2 

R sin rj 2 


(98) 
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COS 02 


X) 

R sin r \2 


(99) 


The relationship between the earth-centered coordinates, X*, X 2 , 

and X 3 , and the perigee system may be found from equations 47 
through 49. Substituting this information along with the fact that 


xp - R cos v 
yp = R sin v 


Zp - 0 


."I 


in the plane 
of orbit 


into (97), (98), and (99) gives 


( 100 ) 


cos r |2 “ 


s in 0 2 


+ 


and 


cos 0 2 - 


(sin i sin u>) cos v + (sin i cos co) 
^sin £2 cos co + cos £2 cos i sin 00 j 
£cos £2 cos i cos co - sin £7 sin 00 j 


|cos £2 cos 00 - sin £2 cos i sin go 
|cos £2 sin go + sin £2 cos i cos co 


sin v 


cos v 
sin ri 2 

sin v 
sin t ) 2 


COS V 
sin r\ 2 

sin v 
sin r \ 2 


( 101 ) 


( 102 ) 


(103) 


The sine of r ) 2 is given by 


sin t \ 2 


Vl - (cos X] Z ) Z 
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since r )2 always lies between zero and n. 

The components of the other vectors needed in the calculations 
of the view factors may be determined from the above equations and 
data. To integrate the view factors the areas A Z j must be chosen. 
This selection is^asy for the cylindrical part of the satellite since 
the direction of N z does not change with axial distance (see equation 
7 3) but only with angular position, £ . Therefore, if the cylinder is 
approximated by a circumscribed polygon, N 2 will have a constant 
direction over the surfaces of any side. The hemispherical ends 
present more of a problem since the unit normal varies with both 
angular coordinates. By approximating the hemispherical surface 
with tangent planes the surface normal is made constant over the 
areas. Using methods such as these, the surfaces of a satellite may 
be approximated to any degree of accuracy desired by choice of the 
sizes of area, A 2 p The resulting integrals have been programmed 
in a general manner on the IBM 7090, using the numerical method of 
Gauss. Typical curves for the planetary, albedo and solar radiation 
view factors are shown in FIG 18 through 20 for each orientation. 

A detailed description of the irradiation may now be determined, 
using equations 13, 17, and 24. A typical curve showing the earth 
planetary radiation incident on the cylindrical part of the satellite, 
at t, = 0, is presented in FIG 21. The solid line gives the heat flux 

for a non- spinning vehicle and the dotted line for a spinning satellite. 
The total incident energy on the cylinder is listed in Table 3 according 
to the source of the radiation and the orientation of the satellite. The 
same information is given in Tables 4 and 5 for the positive and 
negative hemispheres, respectively. 


RESULTS AND DISCUSSION 


A check on the view factor equations, as derived in this report, 
may be obtained for certain simple cases. Smolak in Reference 3 
gives the view factor equations for planetary radiation to a flat plate 
parallel or perpendicular to the local horizon. If the plate is 
perpendicular to the local horizon this equation is 


F v = 3 


tan' 


'( 


pp/h 


yr + 2 PP /h 


(p P /h) yT 


+ 2pp /hi 


(1 + pp/h)^ 


( 105 ) 
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Figure 18. Typical View Factor Curves for Albedo Radiation to Cylinder Oriented Tangent to the Flight 
Path at Perigee 
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Figure 19. Typical View Factor Curves for Solar Radiation to Cylinder Oriented Toward the Earth 
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Figure 20. Typical View Factor Curves for Earth Radiation to Cylinder Oriented toward the Sun 
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TABLE 3 
Cylinder 


Incident Thermal Energy for One Orbit 


_ 3 _ 

TTDL 


Btu 

ft 2 


Source of 
Radiation 


Ear th 


Tangent Flight 
Path at Perigee 

35 


Sun 


123 


Albedo 


21 


Satellite Orientation 


Toward Earth 


32 


111 


19 


Toward Sun 


36 


0 


22 


Total 


179 


162 


58 
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TABLE 4 

Positive Hemisphere 


Source of 
Radiation 

Incident Thermal Energy for One Orbit — 

27TR^ ft? 

Satellite Orientation 

Tangent Flight 
Path at Perigee 

Toward Earth 

Toward Sun 

Earth 

36 

58 

36 

Sun 

112 

64 

191 

Albedo 

20 

35 

16 

Total 

168 

157 

243 
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TABLE 5 


Negative Hemisphere 


Source of 
Radiation 

Incident Thermal Energy for One Orbit Btu 

2TTR 2 ft 2 

Satellite Orientation 

rangent Flight 
Path at Perigee 

Toward Earth 

Toward Sun 

Earth 

24 

14 

36 

Sun 

86 

131 

0 

Albedo 

24 

8 

27 

Total 

134 

153 

63 
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(The horizontal view factor has previously been given as equation 27.) 
In FIG 22 the results of the numerical integration of the earth view 
factor are shown for a flat plate oriented tangent to the flight path 
at perigee. The values given by the Smolak's formulas are shown as 
circle symbols in FIG 22 and agree exactly with the predictions of 
this report. 

To determine the intensity of the earth's planetary radiation, a 
uniform black body temperature is assumed for the entire surface. 

This assumption is necessary since the local temperature is a function 
of many variables and is not predictable. In Reference. 8, Francis 
made a more detailed analysis of the earth T s radiosity. Using the 
mete r ological data of other researchers, a mean environment was 
obtained for each locale, at various times of the year, and the 
resulting irradiation of a plane surface was studied. It was found 
that the local variations from the uniform temperature results were 
quite high for polar orbits (-17% to +33%) but agreed within 5% for 
equatorial orbits and average conditions on the earth. 

The benefits to be gained from controlling the incident energy 
through vehicle orientation may be determined by studying the results 
of the example problem. By orienting the satellite toward the sun, 
the total heat flux to the cylinder is reduced by 120 Btu/ft 2 each orbit. 
However, the incident energy to the positive hemisphere is increased 
by 76 Btu/ft 2 and the heat load to the negative hemisphere is reduced 
by 50%. Consequently, it is very advantageous from a thermal view- 
to orient the satellite toward the sun. Orienting the vehicle 
toward the earth does not offer much, since the proximity of the earth 
reduces the effectiveness of the satellite as a shield. 

The methods outlined in this report may be extended easily to 
predict the incident energy on satellites of other celestial bodies. The 
major problems would be locating the orbital plane and sun relative 
to the body. However, the extension to a lunar satellite is very simple 
since the sun and the vernal equinox have almost exactly the same 
direction at the moon as at the earth. 
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APPENDIX A 
VIEW FACTORS 


The standard equation for the radiation emitted by black body Aj 
which falls on the body A 2 is 


Qi 


<r T 1 


J J 


cos Yi cos Y 2 


l-Tz 


dA, dA? 


( 1 A) 


The integral in this equation is usually combined with either the 
area of the receiver, or the emitter, and is referred to as the view 
factor. If the area of the emitter is used, the physical meaning of the 
view factor is more obvious. Multiplying the numerator and denominator 
of the R.H. S. of the equation by Aj gives 


Q1-2 = ^T , 4 Aj F lz 


(2A) 


where 


F 12 



cos Yi cos Y 2 



dAj dA 2 


(3A) 


The term 0"Ti 4 Aj is the total energy radiated by area Aj and, 
therefore, the view factor represents the fraction that is intercepted 
by A 2 . However, it is equally correct to base the view factor on the 
area of the receiver, if the numerator and denominator of the R.H.S. 
of equation 1A are multiplyed by A 2 . Then 


o' Ti 4 A 2 F 21 


Ql-2 


(4A) 
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where 


21 


1 r 

f 

cos Yi cos Y 2 

7T A 2 J 

Ai 

J 

A z 

Kzl* 


dA, dA 2 


(5A) 


and the order of the subscripts on the view factor is reversed to 
indicate that the reference used in the definition is the area A 2 . A 
comparison of equations 5A and 3A shows the following relationship 
to exist: 


A 2 F 21 - Aj F )2 (6A) 

This is known as the reciprocity theorem. 

In this report the area of the receiver has been used to define 
the view factor and it is assumed that the integration is piecewise 
independent of the area A 2 . That assumption simplifies equation 5A 
to 


F 21 



cos Yi cos Y 2 



dA x 


(7A) 


The equation for the thermal radiation received is unchanged, retaining 
the form of equation 4A. 

Another identifying property of the view factor is its numerical 
value which ranges from zero to plus one. It should also be recognized 
from studying the form of the radiation heat transfer equation, since 
the product of the intensity of the source and the view factor gives the 
incident radiant energy per unit area and time. 
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APPENDIX B 

COORDINATE TRANSFORMATION 


Using the transformation equations 47 through 49, the unit vectors 
of the satellite triad may be expressed in terms of the geocentric co- 
ordinate system. The method consists of successively applying these 
equations until the desired result is obtained. Consider the ip vector. 


ip - cos to (i 10 ) + sin to (i u ) + (0) (i 12 ) 


or substituting for i 10 , and i n , then 


ip ~ c 


os to [(1) (f 7 ) + (0) (J) + (0) (5)] + 
lin to ^(0) (i 7 ) + cos i (i 8 ) + sin i (i 9 )J 


(IB) 


(2B) 


Expressing i 7 , i 8 , i 9 in terms of earth triad gives 


in — C O S CO 


sin 


J^cos £2 (ii) + sin £2 (i 2 ) + (0) (i 3 )| + 

co cos i | - sin £2 (ix) + cos £2 (i 2 ) + (0) (i 3 )j + (3B) 

sin co sin i £(0) (i x ) + (0) (i 2 ) + (1) {i 3 )j 

which simplifies to 

ip = | cos co cos £2 - sin £2 sin co cos ij (ix) + 

J cos co sin £2 4 sin co cos i cos £2 J (i 2 ) + (4B) 

Jsin co sin ij (i 3 ) 
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The remaining unit vectors are derived similarly and are obtained 


as 


Jp 


sin oo cos £2 - cos oo cos i sin £2 J (ii) + 
^cos oo cos i cos £2 - sin oo sin £2 j (^) + 

|cos oo sin ij (i 3 ) 


and 


j^sin i sin £2 j (ii) + |- sin i cos £2 J ( 12 ) + 

| cos ij (i 3 ) 


(SB) 


(6B) 


The earth-oriented triad may be expressed in terms of the perigee 
system by similar methods. The results are 


O. 

= 

- sin V (ip) + cos v (jp) 

(7B) 

Je. 0 . 

= 

- cos v (ip) - sin v (j p ) 

(8B) 

k E. O. 

= 

k P 

(9B) 

The equations for the sun-oriented triad are 


i s. 0 . 

= 

- sin r (ip) - cos T (j p ) 

(10B) 

Js. 0 . 

= 

sin £ cos r (ip) - sin £ sin T (jp) + cos 

e (k p ) 




(1 IB) 

k s. 0 . 

. 

— ► ^ 

- cos £ cos F (ip) + cos £ sin T (jp) + sin 

i (kp) 


( 12B) 
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APPENDIX C 

SATELLITE SURFACE UNIT NORMALS 


The surface unit normals for a cylindrical satellite with hemi- 
spherical ends, oriented toward the earth, are given below 

Cylinder: 




N* = 

- sin v | 

r 

cos £ 

|cos go cos £2 - sin go cos i sin (ij + 



cos C, 

|cos go sin £2 + sin go cos i cos £2 j (i 2 ) + 



o 

o 

\T* 

sin go sin i J (i 3 )J + 


COS V | 

| cos t, 

sin go cos £2 - cos go cos i sin £2 J (ij) + 



COS £> 

|cos go cos i cos £2 - sin go sin £2j (i 2 ) + 



cos C, 

r i -l (1C) 

cos go sin i (i 3 ) f + 


sin £ 

^ sin i sin £2 | (i 2 ) + sin £ sin i cos £2 J (i 2 ) 


sin £ £ cos i j (i 3 ) 

Positive Hemisphere: 
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-sinvsin|isinp[[cosGocos£2 - sin go cos i sin £2 j ( i i) + 

[ cos co sin £2 + sin go cos i cos £2 J (i 2 ) + 
[ sin go sin i J (i 3 ) J + 

- cos v cos p [ [cos go cos £2 - sin go cos i sin £2 J (i^ + 

[cos go sin £2 + sin go cos i cos £2 J (i 2 ) + 

[sin go sin i J (i 3 )j + (2C) 

- sin v cos p [ [- sin go cos £2 - cos go cos i sin £2 J (ii) + 

|cos go cos i cos £2 - sin go sin £2 j (i 2 ) + 

|cos go sin i J ( is) J + 

cos p. sin p | ^ sin i sin £2 1 (ii) + [ - sin i cos £2 j (i 2 ) + 

[cos i ] (i 3 ) j 


Negative Hemisphere: 


N 2 = 


cos v sin 4* sin <J> [ [- sin go cos £2 - cos go cos i sin £}j (ii) + 

[cos go cos i cos £2 - sin go sin £2j (i 2 ) + 
[cos go sin ij (is) J + 

- sin v sin vJj sin 4> [ [cos go cos £2 - sin go cos i sin £2 j (ii) + 

[cos go sin £2 + sin go cos i cos £2 j (i 2 ) + 
[sin go sin i j (i 3 ) j + 
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- cos v cos <j) | |cos oo cos 12 - sin go cos i sin 12 j (ij) 
[cos oo sin 12 + sin oo cos i cos 12 j (i 2 ) 


[sin co sin i j ( i 3 ) J + (3 

- sin v cos <£> [[-sin oo cos 12 - cos oo cos i sin (i\) 

[cos oo cos i cos 12 - sin oo sin 12 j (i 2 ) 

[cos oo sin i ] u7) ) + 

cos ^ sin cf> | [sin i sin 12 j (ij) + [- sin i cos 12 J (i 2 ) 

[cos i] (i 3 ) j 

The equations for the sun-oriented vehicle are given below. 


- sin F cos C, [ [cos oo cos 12 - sin oo cos i sin 12 J (ij) + 

[cos oo sin 12 + sin oo cos i cos 12 j (i 2 ) + 

[sin oo sin i j (i 3 ) J + 

- cos r cos £, [[- sin oo cos 12 - cos oo cos i sin 12 j (ij) + 

[cos oo cos i cos 12 - sin oo sin 12 ] (i 2 ) + 

[cos 00 sin i j (i 3 ) J + 
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- sin £ cos £ cos F |[cos to cos £2 - sin to cos i sin £2 

|cos to sin £2 + sin <o cos i cos £} 
[sin co sin ij (i 3 ) J + 

sin £ cos £ sin F | J - sin co cos £2 - cos co cos i sin S 

[cos co cos i cos £2 - sin co sin £2j 
[ cos co sin i j (i 3 ) J + 

sin £ sin £ j[sinisin£2j(ii) + | - sin i cos £2 J (i 2 ) 

[cos i j (i 3 ) J + 

Positive Hemisphere: 


N 2 - 


- sin F sin jjl sin p | [cos co cos £2 - sin co cos i sin £2j 

[cos co sin £2 + sin co cos i cos £2 

[sin to sin i j (i 3 ) J + 

- cos r sin |jl sin p [ [ - sin co cos £2 - cos co cos i sin S 

Jcos co cos i cos £2 - sin co sin ftj 

[cos co sin i^ (i 3 ) J + 

sin £ cos F cos p | [cos co cos £2 - sin co cos i sin £2 

[cos co sin £2 “t sin co cos i cos £2 

[ sin co sin i j (i 3 ) J + 


(ii) + 

] (3 + 

(4C) 

(ii) + 
(12) + 


(ii) + 

] <iz) + 

^ (3 + 
( 5 ) + 

| (3 + 

] <3 + 
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cos | cos p [ [ sin i sin £2 j (i 3 ) + [- sin i cos ftj (i 2 ) + 

[cos i j (i 3 ) J + (5C) 

- sin £ sin F cos p | J - sin go cos ft - cos go cos i sin ft j (i x ) + 

Jcos go cos i cos ft - sin go sin ftj (i 2 ) + 

| cos go sin ij (i 3 ) J + 

- cos £ cos r cos [j. sin p | jcos go cos ft- sin go cos i sin ftj (ij) + 

Jcos go sin ft + sin go cos i cos ftj (i 2 ) + 
Jsin go sin i j (is) j + 

cos i sin F cos \x sin p j sin go cos ft - cos go cos i sin ftj (i x ) + 

Jcos go cos i cos ft - sin u sin ftj (i 2 ) + 
£cos go sin i ](S>) + 

sin £, cos (x sin p [ [sin i sin ftj (ij) + [ - sin i cos ftj (i 2 ) + 

[cos i j (i 3 ) J + 

Negative Hemisphere: 


- - sin 


F sin 4 J sin cj> | Jcos go cos ft - sin go cos i sin ftj ( i x ) 
Jcos go sin ft + sin go cos i cos ftj (i 2 ) 
Jsin go sin i j (i 3 ) J + 
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- cos 


[ 


COS 00 COS 


r sin 4* sin <{>|[-sln oo cos ft - cos oo cos i sin ftj (ij) + 

i cos ft - sin oo sinftj (i 2 ) + 
[cos oo sin i j (is) J + 

sin £, cos r cos <j> [[cos oo cos ft - sin co cos i sin ftj (i\) + 


[cos oo si 


in ft + sin oo cos i cos ft 


] (3 


cos 



[sin oo sin i j 

(S>] + 


♦ 11 

sin i sin ft j (ii) 

+ sin i cos 

<3> + 


[cos i j (i 3 ) j + 


(6C) 


- sin £ sin T cos <j> | [- sin oo cos ft - cos oo cos i sin ftj (ij) + 

[cos oo cos i cos ft - sin oo sin ftj (i 2 ) + 

[cos oo sin ij (is) J + 

- cos ^ cos r cos ^ sine}) [ [cos oo cos ft - sin oo cos i sin ftj (ij) + 

[cos oo sin ft + sin oo cos i cos ftj (i 2 ) + 

[sin oo sin i j (i 3 ) J + 

cos £ sin P cos vJj sincj> | [- sin oo cos ft - cos oo cos i sin ftj (ij) + 

[cos oo cos i cos ft - sin oo sin ftj (i 2 ) + 

[cos oo sin ij (i 3 ) J + 
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sin £ 

, cos ^ sin cf> 

| [sin i sin S2 J 

(ii) + 

£ - sin i cos £2 j 



[cos i] (i 3 ) j 
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